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Abstract 

Using 13.5 fb^ 1 of e + e~ annihilation data collected with the CLEO II detector, we have observed 
a narrow resonance decaying to D* + tt°, with a mass near 2.46 GeV/c 2 . The search for such a state 
was motivated by the recent discovery by the BaBar Collaboration of a narrow state at 2.32 GeV/c 2 , 
the L>*j(2317)+, that decays to D+tt°. Reconstructing the D+tt° and D* + tt° final states in CLEO 
data, we observe peaks in both of the corresponding reconstructed mass difference distributions, 
AM(D s tt°) = M(D s tt°) - M(D S ) and AM(D*vr°) = M(D*ir°) - M(D*), both of them at values 
near 350 MeV/c 2 . We interpret these peaks as signatures of two distinct states, the D*j(2317) + 
plus a new state, designated as the D s j (2463) + . Because of the similar AM values, each of these 
states represents a source of background for the other if photons are lost, ignored or added. A 
quantitative accounting of these reflections confirms that both states exist. We have measured the 
mean mass differences (AM(D s tt )} = 350.0 ± 1.2 [stat.] ± 1.0 [syst.] MeV/c 2 for the D*j(2317)+ 
state, and (AM (D*ir )} = 351.2 ± 1.7 [stat.] ± 1.0 [syst.] MeV/c 2 for the new L> sJ (2463) + state. 
We have also searched, but find no evidence, for decays of the two states via the channels D* + j, 
Df'y, and Dfir + TT~. The observations of the two states at 2.32 and 2.46 GeV/c 2 , in the Dfir 
and D* + ir° decay channels respectively, are consistent with their interpretations as cs mesons with 
orbital angular momentum L = 1, and spin-parity J p = + and 1 + . 

PACS numbers: 14.40.Lb, 13.25.Ft, 12.40.Yx 



*Submitted to Physical Review D; supersedes hep-ex/0305017. 
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I. INTRODUCTION 



The BaBar Collaboration has recently reported [1] evidence for a new narrow resonance 
with a mass near 2.32 GeV/c 2 , which decays to D+ir°. The BaBar data are consistent with 
the identification of this state as one of the four lowest-lying members of the cs system 
with orbital angular momentum L — 1, and provisionally it has been named the _D*j(2317) 
meson. A natural candidate would be the 3 P cs meson with spin-parity J p = + , but other 
possibilities, including exotic states, are not ruled out. In this paper, we report on a search 
for the _D*j (23 17) meson, as well as other, possibly related states, in data collected with 
the CLEO II detector in symmetric e + e~ collisions at the Cornell Electron Storage Ring, at 
center-of-mass energies ^/s ~ 10.6 GeV. 

The spectroscopy of P-wave cs mesons is summarized in Ref. [2]. Prior to the observation 
of the _D*j(2317), theoretical expectations [3-7] were that: (1) all four states with L = 1 
are massive enough that their dominant strong decays would be to the isospin-conserving 
DK and/or D*K final states, (2) the singlet and triplet J p = 1 + states could mix, and (3) 
in the heavy quark limit, the two states with j = 3/2 would be narrow while the two with 
j = 1/2 would be broad, where j is the sum of the strange quark spin and the orbital angular 
momentum. Existing experimental evidence [8, 9] for the narrow _D s i(2536) and _D*j(2573) 
mesons which decay dominantly to D*K and DK respectively, and the compatibility of the 
D*j (2573) with the J p assignment as 2 + support this picture. 

The observation by BaBar [1] of the new state at a mass of 2.32 GeV is surprising because: 
(1) it is narrow (with intrinsic width T < 10 MeV), (2) it has been observed in the isospin- 
violating D s tt channel, and (3) its mass (2316.8 ± 0.4 [stat.] MeV/c 2 ) is smaller than most 
theoretical predictions for a + cs state that could decay via this channel. However, points 
(1) and (2) would be obvious consequences of the low mass, since the D^K decay modes 
are not allowed kinematically. We also note that at least two theoretical calculations [10, 11] 
prior to the P*j(2317) + observation had suggested that, in the heavy quark limit, the j = 1/2 
states with J p = + and 1 + could be thought of as chiral partners of the D s and D* mesons, 
and thus would be relatively light. In one model [11] it was proposed that the mass splittings 
between the + and 0~ states of heavy flavored mesons could be as small as 338 MeV/c 2 , 
which is near the _D*j(2317) + — D+ mass splitting of 348.3 MeV/c 2 measured by BaBar. 

Since the initial observation, a number of explanations have appeared [12-19]. Cahn 
and Jackson [12] apply non-relativistic vector and scalar exchange forces to the constituent 
quarks. Barnes, Close and Lipkin [13] consider a quark model explanation unlikely and 
propose a DK molecular state. Similarly, Szczepaniak [16] suggests a Dn atom. Also going 
beyond a simple quark model description, Van Beveren and Rupp [14] present arguments 
for a low mass + cs state based on a unitarized meson model, by analogy with members 
of the light scalar meson nonet. Bali [19] reports on lattice QCD calculations that predict 
signficantly larger + — 0~ meson mass splittings than what has been observed for the 
£^(2317) - D s splitting. 

On the contrary, Bardeen, Eichten and Hill [15] couple chiral perturbation theory with a 
quark model representation in heavy quark effective theory, building on the model described 
in Ref. [11]. They infer that the P*j(2317) is indeed the + cs state expected in the quark 
model, predict the existence of the 1 + partner of this state with a 1 + - T mass splitting 
equal to the + — 0~ mass splitting, and compute the partial widths for decays to allowed 
final states. Godfrey [17] and Colangelo and De Fazio [18] find that the radiative transistion 
of the P*j(2317) should be significant if it is indeed a cs state. 
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The goals of the analysis presented here are to use CLEO data to provide independent 
evidence regarding the existence of the _D*j (231 7), to shed additional light on its properties, 
and to search for decays of other new, possibly related states. In particular, we address 
the following questions. Are the electromagnetic decays D s j or _D*7 observable in light of 
the isospin suppression of the strong decay to D s tt° ? Are other strong decays observable 
such as D*7T°, or the isospin-conserving but Okubo-Zweig-Iizuka (OZI) suppressed [20] decay 
D s -n + it~ ? If the -D*j(2317) is the expected + cs state, might the remaining 1 + state also 
be below threshold for decay to D*K, as suggested in Ref. [15], and thus be narrow enough 
to be observable in its decays to D*ir°, D s ^ or -D*7? 

This article [21] is organized as follows. After describing the detector and data set in 
Section II, we summarize the reconstruction of the _D*j(2317) + — > D+tt° decay channel 
in Section III, including efforts to understand and exclude contributions from known back- 
ground processes. We then report in Section IV on searches for other possible decay channels 
as described in the preceding paragraph. In Section V, we report on the appearance of a 
statistically significant signal in the D* + n° channel at a mass of 2.463 GeV/c 2 , not compat- 
ible with a kinematic reflection of the _D*j(2317) + . We describe a quantitative analysis of 
the signals in the Dfir and D* + ir° channels, leading us to infer the existence of two distinct 
states. Based on this conclusion, we discuss the properties of these two states in Section VI, 
after which we summarize the principal results of the analysis. 

II. DETECTOR AND DATA SET 

The analysis described here is based on 13.5 fb _1 of e + e _ collision data collected between 
1990 and 1998. CLEO II is a general purpose, large solid angle, cylindrical detector featuring 
precision charged particle tracking and electromagnetic calorimetry, and is described in detail 
in Refs. [22, 23]. In its initial configuration, the tracking system was comprised of a six-layer 
straw tube chamber just outside of a 3.2 cm radius beryllium beam pipe, followed by a 10 
layer hexagonal cell drift chamber and a 51 layer square cell drift chamber, immersed in a 
1.5 T magnetic field generated by a superconducting solenoid. In 1995, the beam pipe and 
straw tubes were replaced by a 2.0 cm radius beam pipe plus three layers of silicon strip 
detectors each with double-sided readout, and a helium-propane gas mixture replaced the 
argon-ethane mixture previously used in the main drift chamber. 

Beyond the tracking system, but within the solenoid, were also located a 5 cm thick plas- 
tic scintillation counter system for time-of-flight measurement and triggering, and a barrel 
calorimeter consisting of 6144 tapered CsI(Tl) crystals 30 cm in length, arrayed in a pro- 
jective geometry, with their long axis oriented radially with respect to the e + e~ interaction 
point. An additional 1656 crystals were deployed in two end caps to complete the solid angle 
coverage. The excellent energy and angular resolution of the calorimeter is critical for the 
reconstruction of 7r° — > 77 decays as well as single low-energy photons such as those emitted 
in the D* + — > D.^ transition. 

III. CONFIRMATION OF D*j(2317)+ D+tv° 

The search for the -D*j(2317) was carried out by reconstructing the D+tt° state, using the 
D+ — > (f)-K + channel with — > K + K~. Charge conjugation is implied throughout this article. 
Pairs of oppositely charged tracks were considered as candidates for the decay products of 



4 



the if the specific ionization (dE/dx) was measured in the main drift chamber to be within 
2.5 standard deviations of the expectation for a kaon, and if the invariant mass of the K + K~ 
system was within ±10 MeV/c 2 of the mass. A third track with dE/dx consistent with 
the expectation for a pion was combined with the K + K~ system to form a candidate 
with mass M{KKn). To improve resolution we adjust the momenta of the three particles 
subject to the constraint that their trajectories intersect at a common point corresponding 
to the decay point of a D s meson. When fitted to a Gaussian, the observed mass peak 
has a standard deviation (a) of 6.5 ± 0.4 MeV/c 2 in our data, consistent with CLEO Monte 
Carlo simulations of D s production and decay plus a GEANT-3 [24] based simulation of 
particle propagation and detector response. 

Clusters of energy deposition in the calorimeter unassociated with charged particle in- 
teractions were identified as potential photon candidates. To be considered as candi- 
dates for the photons from ir° — > 77 decay, clusters with energy greater than 100 MeV 
located in the central region of the calorimeter (|cos#| < 0.71, where 6 is measured 
with respect to the beam axis) were selected. Pairs of photons were required to satisfy 
—3.0 < [M(77) — M 7r o]/cr(77) < 2.5, where M(77) is the invariant mass of the two photons 
and 17(77) is the expected resolution on this mass. For each cluster being considered as a 
photon candidate, we additionally required that the lateral profile of energy deposition in the 
calorimeter be consistent, at the 99% confidence level, with expectations for photons. This 
requirement removes spurious photon candidates that are mainly due to inelastic interac- 
tions of charged hadrons or long-lived neutral hadrons. The peak in the M(77) distribution 
for photon-pairs accompanying a D+ candidate with M(D S ) = M(KKir) between 1.9565 
and 1.9805 GeV/c 2 has a = 5.8 ±0.4 MeV/c 2 in our data, consistent with expectations from 
the Monte Carlo simulations. Once identified as a ir° candidate, the directions and energies 
of the two photons are adjusted with a kinematic fit to reconstruct to the known value [8] 
for the 7r° mass M w o. 

To suppress combinatoric backgrounds, we further required that the momentum of the 
-D+7T candidate be greater than 3.5 GeV/c. We also required that the helicity angle 
of the — > K + K~ decay satisfy the requirement |cos#/i| > 0.3, where Oh is the angle 
between the K + momentum vector measured in the rest frame, and the momen- 
tum vector measured in the D s rest frame. The expected distribution from real de- 
cays varies as cos 2 9 h , whereas combinatoric backgrounds tend to be flat. For D s n° com- 
binations satisfying these requirements, we plot the mass M(D s 7r°) = M(KK7iir°) and 
the mass difference AM(D s tt°) = M(D s tv°) - M(D S ) in Fig. 1(a) and (b), respectively. 
To improve the experimental resolution on M(D s tc°), the known value of the D s mass, 
Md s = 1968.5 ± 0.6 MeV/c 2 [8], has been used to determine the energy of the KKir system 
from its measured momentum in Fig. 1(a); this substitution is not done for AM(D s 7r ) in 
Fig. 1(b), or for the calculation of other mass differences entering this analysis. 

The narrow peaks in Fig. 1 at a mass near 2.32 GeV/c 2 and a AM(D s ir°) near 350 
MeV/c 2 are in qualitative agreement with the BaBar observation. We note that there are 
no peaks in this region when KKn combinations with M(KKtt) lying in D s side band 
regions are combined with a 7r°. The other feature of note in the spectra is the sharp signal 
from D* + — > D+7i° [25] near the kinematic threshold. In addition, Monte Carlo simulations 
of inclusive multi-hadron production via e + e~ — > qq (q = u, d, s, c) give M(D s 7r°) and 
AM(D s 7i°) spectra that reproduce the features observed in the data, except for the peaks 
near 2.32 GeV/c 2 and 350 MeV/c 2 in the respective plots. This is also illustrated in Fig. 
1, where the normalization for the qq Monte Carlo spectra is fixed by the ratio of the 
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FIG. 1: Distributions of (a) the masses M(D s tt°) of the D s ir° candidates and (b) the mass differ- 
ences AM(D S 7T°) = M(D s tt°) — M(D S ) for events satisfying cuts on M(KKtt) consistent with the 
D s mass and M^yy) consistent with the 7r° mass, as described in the text. The points represent 
the CLEO data, while the solid histogram is the predicted spectrum from the Monte Carlo simu- 
lation of e + e~ — ► qq events. The predicted spectrum is normalized absolutely by the ratio of the 
equivalent luminosity of the Monte Carlo sample used to the luminosity of the CLEO data sample. 
The overlaid curve represents the results from a fit of the data to a Gaussian signal function plus 
a second-order polynomial background function. 

luminosity of the data sample to the equivalent luminosity of the Monte Carlo sample. This 
normalization is known to a precision of approximately ±5%. 

The agreement between the Monte Carlo and data distributions in Fig. 1 in normalization 
as well as shape demonstrates that the simulation of 'random' photons accompanying D s 
decays is accurate. The accuracy of this simulation is important for our detailed analysis of 
this signal, described in Section VA. 

We have investigated mechanisms by which a peak at 2.32 GeV/c 2 could be generated 
from decays involving known particles, either through the addition, omission or substitution 
of a pion or photon, or through the mis-assignment of particle masses to the observed 
charged particles. In no cases were narrow enhancements in the M(D s 7r°) spectrum near 
2.32 GeV/c 2 observed. We will discuss the issue of backgrounds from a new resonance at 
2.46 GeV/c 2 when we describe our studies of the D* + ir° final state. 
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From a binned maximum likelihood fit of the AM(D s 7r°) distribution to a Gaussian signal 
shape and second-order polynomial background function, we obtain a yield of 165 ±20 events 
in the peak near 350 MeV/c 2 . In this fit, the mean and Gaussian width of the peak are 
allowed to float. These parameters are determined to be {AM(D s ir )) = 349.4±1.0 MeV/c 2 
and a = 8.0t\i MeV/c 2 , where the errors are due to statistics only. The peak is somewhat 
broader than the expected mass resolution of 6.0 ± 0.3 MeV/c 2 , determined from Monte 
Carlo simulations. The detection efficiency associated with the reconstruction of the full 
D* sJ (2317) + -> D+tt , D+ -> <M+, -> K + K- decay chain is (9.73 ± 0.57) % for the 
portion of the _D*j(2317) + momentum spectrum above 3.5 GeV/c, where this efficiency does 
not include the D s and <fi decay branching fractions. 

Thus, we confirm the existence of a peak in the D s ir° mass spectrum that cannot be 
explained as reflections from decays of known particles. Our measurements of the mean 
mass difference and width of the peak are consistent with the values obtained by BaBar [1] 
for the .D*j(2317) + . Further discussions of the width, as well as of systematic errors in the 
measurements of the mass and width of the _D*j (231 7) appear later in this article. 

IV. SEARCHES FOR D* J (2317) IN OTHER CHANNELS 

The conclusion that the _D*j(2317) is a new narrow resonance decaying to D s n° leads 
to two questions: (1) are there other observable decay modes, and (2) might additional 
new cs resonances also exist in which normally suppressed decay modes such as D^n° are 
dominant? To answer these questions we have searched in the channels D s j, D*j, D*n° 
and D s it + it~ . 

If the _D*j (23 17) is a + L = 1 cs meson, as has been suggested [15], it could decay via 
S- or D-wave to D*j, but would not be able to decay to D s ^ due to parity and angular 
momentum conservation. Consequently, observation of one or both of these channels would 
be interesting. On the other hand, if neither channel is seen, this would not be too surprising 
since these are electromagnetic decays, and the D s n° decay, while isospin violating, is not as 
severely phase-space suppressed as in the case of the corresponding decay of the D* where 
the electromagnetic decay dominates. The BaBar data show no evidence for either channel, 
however no upper limits were reported on the branching ratios for these channels. 

With regard to strong decays, the D s n + n~ final state is kinematically allowed and isospin 
conserving, but would be suppressed by the OZI rule. This is in contrast to the D s n° channel 
for which one mechanism would be decay to a D s plus a virtual 77, with production of the 
7T° via i]-7i° mixing [26]. However, angular momentum and parity conservation forbid the 
decay of a + state to three pseudoscalars. Thus, observation of the D s n + 7i~ channel would 
be strong evidence against the interpretation of the _D*j(2317) as a + meson. 

Finally, it is possible that the remaining L = 1 cs state with J p = 1 + could also be 
light enough that decays to D*K would be kinematically forbidden. In this case, the strong 
isospin-violating decay of this 1 + state to D*tt° could occur via S-wave (the electromagnetic 
decays to -D s 7 or D*^ would also be possible), and thus a narrow peak in the AM (£>>°) = 
M(D*ir°) — M(D*) spectrum would be a signature of such a state. 
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FIG. 2: The mass difference AM(D s 7nr) = M(D s Tnr) - M(D S ) for candidates, as de- 

scribed in the text. 



A. Searches for D*j(2317)+ decays to D+-k+7v~, D+-y and D*+-y 

To look for these channels we select events containing D+ — > 07r+ candidates as in the 
D S 7T° analysis. For the D s ii + 7r~ channel, we combine the D s candidates with two oppositely 
charged tracks, and plot the mass difference AM(D s 7m) = M(D s 7T7t) — M(D S ). As shown 
in Fig. 2, no signal is evident in the vicinity of 350 MeV/c 2 . 

To search for states decaying to -D + 7, we have formed -D+7 combinations by selecting 
photons of energy greater than 150 MeV. To select D* + candidates for use in other searches, 
we relax this to include photon candidates with energy above 50 MeV. We ignore photons 
that can be paired with another photon such that M( / -f-f) is consistent with n decay. The 
inclusive AM(D S 7) = M(D S ^) — M(D S ) spectrum for this sample is plotted in Fig. 3(a), 
illustrating that a large D* sample can be obtained. For decay modes with a D* in the final 
state, we select D s '~f combinations where the mass difference AM{D s r y) is reconstructed to 
be between 0.1308 and 0.1568 GeV/c 2 . 

Also visible in Fig. 3(a), are regions of the AM{D s r y) spectrum where decays of the 
_D*j(2317) (or of a possible higher mass state) into D s ^ would appear. There is no evidence 
for a signal near 350 MeV/c 2 corresponding to a M(Dg'y) in the vicinity of 2.32 GeV/c 2 . 

The same conclusion holds for the _D*7 final state, shown in Fig. 3(b), where we combine 
selected D* candidates with photons of energy above 150 MeV. The peak in the AM{D* s r y) 
spectrum in Fig. 3(b) near 150 MeV/c 2 is due to real D* + — > Df^y decays in which a random 
photon has been combined with the candidate to form the D* candidate, and the actual 
photon from this transition is combined with this system to form the D* sJ candidate. There 
is no sign of any structure in this spectrum near 205 MeV/c 2 , where a signal from D* sJ (2?>Yl) 
decay would be expected. 



B. Search for D* aJ (2317)+ decays to D*+ir° 

We have also searched in the D* + ir° channel for D* sJ states. To maintain efficiency for this 
final state, we do not veto D* + candidates where the photon used in the D* + reconstruction 
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FIG. 3: (a) The spectrum of the mass difference AM(D S ~/) = M(D S ~/) - M(D S ), plotted on a 
logarithmic scale. The peak is due to the transition D* + — > Df^y. (b) The spectrum of the mass 
difference AM(L>* 7 ) = M(D*j) - M(D* S ) for D*j candidates. 



can be combined with an extra photon to form a 7r° decay candidate. We also applied slightly 
less restrictive track quality and shower shape criteria than in the D s n° analysis. As with 
the modes involving D* candidates described the preceding section, the energy of photons 
selected for reconstruction of the D* — > D s '-f decay is required to satisfy E 1 > 50 MeV. The 
D*7i° candidates are required to have momenta above 3.5 GeV/c. Fig. 4(a) shows the mass 
difference plot for events with candidate Df — > 07r + , D* + — > decays plus di-photon 
combinations consistent with tc° decay. If the D* sJ (2317) + were to decay to the D* + -k° final 
state, a peak would be expected at a AM(D*ir°) ~ 205 MeV/c 2 . Although we see no 
evidence for such a peak, there is a significant excess in a narrow region near 350 MeV/c 2 . 
We discuss the properties of this new peak in the following section. 



V. OBSERVATION OF A NEW STATE AT 2.463 GeV/c 2 

From a fit to a signal Gaussian signal function plus a polynomial background function, 
we observe a peak in Fig. 4(a), comprised of 55 ± 10 D*7r° combinations, at (AM(D*n )) = 
349.8 ±1.3 MeV/c 2 . The fit yields a Gaussian width of 6.1 ± 1.0 MeV/c 2 for the peak, 
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FIG. 4: (a) The mass difference spectrum AM(D*ir°) = M(D s 'yK ) — M(D S ^) for combinations 
where the D s ~f system is consistent with D* decay, as described in the text, (b) The corresponding 
spectrum where D s 7 combinations are selected from the D* side band regions, defined as 20.8 < 
\AM(D sl ) - 143.9MeV/c 2 | < 33.8MeV/c 2 . 



consistent with our mass resolution of 6.6 ± 0.5 MeV/c 2 . The existence of this peak leads 
us to investigate the possibility of a second narrow resonance with a mass near 2.46 GeV/c 2 
that decays to D* + n°. We note that a similar peak is also present in the M(D*ir°) spectrum 
observed by BaBar [1] , although BaBar does not claim this as evidence for a new state. For 
ease of notation, we refer to the postulated particle as the -D s j (2463) + . 



A. Analysis of cross feed between D+7T and D* + 7T° samples 

The kinematics of the Dfir and D* + 7r° decays are quite similar, and it is possible that 
they can reflect into one another. For example, by ignoring the photon from the D* decay 
in D s j(2A63) + — > D* + tt° decays, nearly all the putative signal combinations form a peak 
in the AM(D s tt°) spectrum in the same region as the D* sJ (2317) signal described in pre- 
vious sections of this article. We refer to the background entering via this scenario in the 
-D*j(2317) + — > D^tt sample as 'feed down' from the _D s j(2463) state. The impact of ne- 
glecting the photon is that this peak in AM(D s n°) is broader than that for real _D*j(2317) 
decays. From Monte Carlo simulations, we determine the width of this smeared peak to be 
a = 14.9 ± 0.4 MeV/c 2 . 
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It is also possible that a Dfn candidate can be combined with a random photon such 
that the -D+ 7 combination accidentally falls in the D* + signal region defined earlier. In this 
case, £>*j(2317) + -> D+ir° decays would reflect or 'feed up' into the D sJ (2463) + -> D*+ir° 
signal region. A Monte Carlo simulation of L>*j(2317) + production and decay to D+ir° 
shows that this does happen, but only for approximately 9% of the reconstructed decays. 
The peak in the AM(D*7i°) distribution generated by this feed up sample is also broadened 
relative to the expectation for real -D s j (2463) decays, analogous to the smearing of the feed 
down kinematics discussed in the preceding paragraph. 

We can extract the number of real -D*j(2317) + — > D+ir° decays reconstructed in our data, 
denoted as Ro, as well as the number of real _D s j (2463) + — > D* + 7i° decays, denoted as R±, 
taking into account that the corresponding real signal decays in one channel can enter the 
candidate sample for the other channel as described above. The following linear equations 
relate the real to observed numbers: 

N = R + f 1 R 1 (1) 

N, =R 1 + f R , (2) 

where N and Ni are the numbers of observed decays in the D s tt° and D*7i° channels 
respectively, and Rq and Ri are the number of real decays produced times the efficiency 
to observe them in the corresponding signal decay channels. The coefficients /o and f\ 
are the feed up and feed down probabilities relative to the reconstruction efficiency for the 
respective signal modes. We note that these relations represent first-order approximations; 
higher-order corrections, such as that due to the scenario where the D s 7i° system from a 
real D s j(2463) decay is combined with an unrelated photon to form a feed up _D<,j (2463) 
candidate, are negligible in the present case. 

The observed number of decays in the D*7r° channel is Ni = 55 ± 10, obtained from the 
fit to the peak in Fig. 4(a) described above. For N , it is desirable to obtain a D s 7r° sample 
selected with criteria that most closely match those used to select D*tt° combinations, and 
that is enriched in _D*j (23 17) decays relative to feed down from _D s j (2463) decays. Thus we 
apply the same selection criteria that were used for the D*n° sample, but without selecting 
the photon from the D* — > D s '~f transition. To measure the event yield in this sample, we 
fit the peak in the AM(D s n ) distribution to a Gaussian with its width fixed to the Monte 
Carlo expectation for _D*j(2317) decays. In this fit, a significant fraction of feed down 
combinations is counted as part of the combinatoric background rather than as signal. We 
obtain N = 190 ± 19 candidates. This sample effectively constitutes the source of potential 
feed up candidates. The difference between this yield and the 165 ± 20 events reported in 
Section III is consistent with the different acceptances for the two sets of selection criteria. 

From the Monte Carlo simulations we measure fo = 0.091 ± 0.007 ± 0.015 for the prob- 
ability that a reconstructed _D*j(2317) — > D s ir° can be combined with a random photon to 
mimic a -D s j (2463) — > D*rr° decay. The first error is due to limited Monte Carlo statistics 
while the second is due to systematic uncertainties associated with: (1) the modeling of extra 
photons in the simulations, and (2) the fraction of such combinations that are counted by 
the fit to the AM(D*ir°) distribution as contributing to the Gaussian signal. As indicated 
above, this fraction counted by the fit is less than one due to the smearing of AM(D*n°) 
that results when an unrelated photon is added to the D s 7r° system. The agreement between 
the data and Monte Carlo distributions in Fig. 1 lends confidence in the modeling of extra 
photons. We assign a relative systematic uncertainty of 5% based on this and on studies of 
combinations entering D* side bands described in the following section. To study the second 
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source of systematic uncertainty, we have carried out fits to the AM(D*7r°) distribution in 
which the width of the Gaussian signal function was fixed to ±1<t relative to the central 
value obtained from the nominal fit. Based on the resulting variation in event yields, we 
have estimated a relative uncertainty on f of 16% from this source. 

We also obtain j\ = 0.84 ± 0.04 ± 0.10 from Monte Carlo simulations, where the first 
error is statistical and the second is due to systematic errors. This includes the probability 
of feed down as well as the photon finding efficiency If all -D s j (2463) — > D*n° decays 
with a reconstructed D s plus 7r° combination were to be counted as D*j(2317) decays, fi 
would simply be one divided by the acceptance for finding the photon from the D* — > D s j 
transition. However, because the AM(D s 7r°) distribution for the feed down background 
is broadened, a significant fraction of these combinations are not counted as part of the 
Gaussian signal, instead being absorbed into the polynomial background. The contributions 
to the relative systematic error on f\ are estimated to be 5% from the uncertainty on the 
photon-finding efficiency and 11% from the uncertainty on the probability of feed down, 
obtained by performing alternate fits to the AM(D s 7r°) distribution. 

Inverting Equations 1 and 2, we find that R = 155 ± 23 decays and R\ = 41 ± 12 decays, 
where the uncertainties include both statistical and systematic sources. The result for i?i 
demonstrates the existence of a state at 2463 MeV/c 2 . The significance of the signal for this 
state, accounting for statistical and systematic errors, is determined to be in excess of 5a 
by computing the probability for the combinatoric background plus the feed up background 
to fluctuate up to give the observed yield in the signal region in Fig. 4(a). 

B. Further evidence for the Z? s j(2463)+ — *■ D*ir° decay 

We conclude from the analysis described in the preceding section that a new state, the 
-D<,j(2463), exists in addition to the _D*j (23 17) state reported by BaBar, because feed up from 
the -D*j(2317) is only a minor background component (~ 25%) of the narrow peak observed 
in Fig. 4(a). To provide further support for this conclusion, we have directly measured 
the feed up background in Fig. 4(a) due to D* sJ (2317) + — > D+ir° plus random photon 
combinations, by selecting combinations in D* side band regions in the D s / ~fTT° sample. 
The M(D s / jtt ) — M(Ds'-f) distribution for this sample, plotted in Fig. 4(b), shows only a 
small enhancement in the region of the _D s j (2463), demonstrating that the background from 
_D*j(2317) decays indeed constitutes only a small fraction of the entries in the _D s j(2463) 
peak. 

We performed a binned likelihood fit of the spectrum in Fig. 4(a) to a Gaussian signal 
shape plus a second-order polynomial plus the spectrum from the D* side band region in 
Fig. 4(b) with its normalization fixed. From this fit, we obtain Ri = 45.7 ± 11.6 decays, 
consistent with the value of Ri obtained from Equations 1 and 2. From the change in the 
likelihood of fits performed with and without the _D s j (2463) signal contribution, we infer 
that the statistical significance of the signal is 5.7 o. 

Finally we note that the width of the peak in Fig. 4(a), a = 6.1 ± 1.0 MeV/c 2 , is consistent 
with the detector resolution. If the origin of this peak was feed up from D*j(2317) + — > D+ir° 
decays, then the effect of including unrelated photons to form D*n° candidates would be to 
smear out the AM(D*n°) distribution, in the same way that the feed down background to 
the _D*j (23 17) state is broadened as described in the preceding section. From fits to Monte 
Carlo simulations of this feed up process, the expectation for the width is determined to be 
a = 14.9 ± 0.6 MeV/c 2 . Thus, the narrowness of the peak in Fig. 4(a) also rules out the 
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possibility that the peak is dominantly due to feed up from _D*j(2317) + decays. 

VI. PROPERTIES OF THE £)* J (2317)+ AND £> sJ (2463)+ STATES 

A. Mass and width of the D* J (2317)+ 

Having obtained evidence for the -D s j (2463) state, and having characterized the back- 
ground that it contributes in the AM(D s ir°) mass difference spectrum, we are now able to 
further address properties of the _D*j(2317) state. We recall that our measurement of the 
width of the peak in Fig. 1 is a — 8-0~L\'.i MeV/c 2 , somewhat larger than our mass difference 
resolution, a = 6.0 ± 0.3 MeV/c 2 . This difference is consistent with predictions from Monte 
Carlo simulations where we include both D s j(2A63) and -D*j(2317) production, since roughly 
18% of the observed Dfn decays in the _D*j(2317) signal region enter as feed down from the 
-D<,j(2463) state, this 'background' peak having an expected width of a = 14.9 ±0.4 MeV/c 2 . 

To better determine the mass and natural width of the _D*j (23 17), we carry out a binned 
likelihood fit of the peak in the AM(D s ir°) spectrum in Fig. 1(b) to a sum of two Gaussians, 
one for the _D*j (23 17) signal, and one to account for the feed down from the _D s j (2463). 
Allowing the means and widths of both Gaussians to float, we obtain (AM(D s -k )) = 350. 0± 

1.2 MeV/c 2 with a = 6.0 ± 1.2 MeV/c 2 for the D* sJ (2317) component. The mean mass 
difference and width for the feed down component are 344.9 ±6.1 MeV/c 2 and 16.5 ± 

6.3 MeV/c 2 , respectively. The errors in the above values are due to statistics only; systematic 
errors are discussed below. Both widths are consistent with predictions from Monte Carlo 
simulations in which the two states are modeled with a natural width of zero. 

We have also carried out fits in which one or both of the widths of the Gaussians were 
fixed to values determined by the Monte Carlo. In all cases the results were consistent with 
the results from the fit described above. We have also tried to obtain a purer D*j(2317) 
sample by vetoing combinations with photons that can be combined with the D s candidate 
to form a D*, thereby removing some of the feed down background from the _D s j (2463). This 
veto marginally improves the D s 7i° signal when we fit with two Gaussians, and the mass 
and width change by only a small fraction of the statistical uncertainty. The systematic 
uncertainty for (AM(D s tt )) receives contributions from uncertainties in the characterization 
of the .D s j(2463) feed down and from uncertainties in the modeling of the energy resolution 
of the calorimeter. We estimate the total systematic error on the mass difference to be 1.0 
MeV/c 2 . Based on these studies, we limit the natural width of the _D*j(2317) to be T < 7 
MeV at the 90% confidence level (C.L.). 

B. Mass and width of the D sJ (2463)+ 

From the fit to the distribution resulting from the subtraction of Fig. 4(b) from Fig. 4(a) 
reported in Section VB, we obtain (AM(D*n )) = 351.2 ± 1.7 ± 1.0 MeV/c 2 for the mass 
difference between the _D s j(2463) and the D*. The first error is statistical and the second 
is the systematic uncertainty which is the same as that presented in the previous section 
for the .D*j(2317) — D s mass difference. From our fits to data and Monte Carlo AM(D*tt°) 
distributions, we also infer a 90% C.L. upper limit on the natural width (r) of the _D s j(2463) + 
state to be 7 MeV. 
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TABLE I: The 90% C.L. upper limits on the ratio of branching fractions for D*j(2317) to the 
channels shown relative to the D+ir° state. Also shown are the theoretical expectations from 
Ref. [15], under the assumption that the D* sJ (2317) is the lowest-lying + cs meson. 



Final State Yield Efficiency Ratio (90% C.L.) Prediction 

-D+7T 135 ± 23 (9.7 ±0.6)% — 

D+7 -19 ±13 (18.5 ±0.1)% < 0.052 

D*+7 -6.5 ±5.2 (7.0 ±0.5)% < 0.059 0.08 

D+TT+ir- 2.0 ±2.3 (19.8 ± 0.8) % < 0.019 

D*+ir° -1.7 ±3.9 (3.6 ±0.3)% < 0.11 



C. Production Properties 

We now give a measure of the production rates of -D*j(2317) and D s j(2463) mesons. A 
full understanding would require the determination of the fragmentation functions of both 
particles and their branching ratios into the final states we observe. To minimize systematic 
errors, we report the relative yields with respect to Df production, where all putative 
charmed-antistrange systems have momenta greater than 3.5 GeV/c. We use all observed 
events for each channel, which includes direct production and any contributions from decays 
of higher mass objects. Then 

a-B(D* sJ (23l7)^Dfn°) 
* (Df) 

a-B(D sJ (2AQ3)^ D* s + n°) 
° (Df ) 

We also note that 

a ■ B(D*+(2112) ->D+7) 

V s K ._J 2-^ = 0.59 ±0.03 ±0.01, (5) 

here and above, the first error includes the statistical and systematic errors on the event 
yields while the second includes the systematic errors for photon detection (2%), and for 7r° 
detection (5%). 

D. Decays of £^(2317) to other final states 

With regard to the alternate D*j(2317) decay channels described earlier, in which no 
signals were observed, we summarize the limits on the branching fractions relative to the 
-D+7T mode in Table I. The normalization for these limits is based on the determination 
that (81.7 ± 5.7) % of the observed yield of 165 ± 20 entries in the peak of the AM(D s n°) 
spectrum in Fig. 1(b) are attributable to D*j(2317) — > D s ir° decay after accounting for the 
feed down from decays of the D sJ (2AQ3) state to D*ir°. We have estimated the systematic 
error on this yield to be ±16 entries by varying selection criteria and the parameterization of 
signal and background shapes used in the fit to Fig. 1. The event yields for the various final 
states are obtained by fitting the mass difference distributions to Gaussians with each mean 
fixed to the result from the Dfn° channel and each width given by the resolution determined 



= (7.9 ± 1.2 ±0.4) x 10~ 2 , 
= (3.5 ±0.9 ±0.2) x 10~ 2 , 



(3) 
(4) 
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TABLE II: The 90% C.L. upper limits on the ratio of branching fractions for D s j(2463) to the 
channels shown relative to the D* + ir° state. Also shown are the theoretical expectations from 
Ref. [15], under the assumption that the D s j(24Q3) is the lowest-lying 1 + cs meson. 



Final State Yield Efficiency Ratio (90% C.L.) Prediction 

D*+tt° 41 ± 12 (6.0 ± 0.2) % — 

D+7 40 ± 17 (19.8 ± 0.4) % < 0.49 0.24 

£>*+7 -5.1 ±7.7 (9.1 ±0.3)% < 0.16 0.22 

D+TT+ir- 2.5 ±5.4 (19.5 ±1.5)% < 0.08 0.20 

D* J (2317)+7 3.6 ±3.0 (2.0 ±0.1)% < 0.58 0.13 



from the simulation of the corresponding decay mode. Uncertainties are dominated by the 
statistical error on the fitted yields and limits on the relative rates are calculated assuming 
a Gaussian distribution with negative values not allowed. 

E. Decays of D s j(2463) to other final states 

Unlike the case of a + state, the D s 7i + n~ decay mode, as well as both radiative decay 
modes D s ^y and D*^ are allowed for a state with J p = 1 + . From fits to the mass difference 
distributions displayed in Figs. 2 and 3 for peaks in the regions where a contribution from 
the -D s j(2463) would appear, we find no evidence of decays to any of these final states. We 
summarize the limits obtained on these decays, relative to D*7r°, in Table II. 

Despite a high relative efficiency, the limit on the decay -D s j(2463) + — > -0^7 is less 
stringent than those on the decays to D* +/ y and D+7t + tt~. This is due to an excess of 
combinations in the signal region. From fits performed with and without the signal Gaussian, 
we determine that the statistical significance of this excess is 2.4 standard deviations. 

If the _D s j(2463) + is a 1 + state, then it is also possible for it to undergo a P-wave 
radiative decay to _D*j(2317) + 7 [27]. We have looked for this transition in our D s ^tt sample. 
To reduce backgrounds from £> s j(2463) + — > -D* + 7, we required that the D s ir° system be 
consistent with the decay of the D* 8j (2ZYI), namely that \AM(D s tt°) - 350.0 MeV/c 2 | < 
13.4MeV/c 2 (~ 2a based on Monte Carlo simulations). We also required that the D s , ~f 
system be inconsistent with D* decay at the la level (the corresponding AM(D S 7) must 
deviate from the expected value for this decay by more than 4.4 MeV/c 2 ), and that the 
momentum of the 7r° be inconsistent with the -D s j (2463) — > D*7r° transistion, also at the la 
level. The M(D s 7r 0/ -f) — M(D s ir°) distribution, plotted in Fig. 5, provides no evidence for a 
signal in the vicinity of 150 MeV/c 2 . Due to the tightness of these cuts, the efficiency for 
detecting this decay is roughly a factor of three smaller than for the D*ir° decay mode. The 
90% C.L. upper limit for this channel is reported in the bottom row of Table II. 

VII. CONCLUSIONS AND DISCUSSION 

In summary, data from the CLEO II detector have provided confirming evidence for the 
existence of a new narrow resonance decaying to Dfir , with a mass near 2.32 GeV/c 2 . 
This state is consistent with being the + member of the lowest-lying P-wave cs multiplet. 
As summarized in Table I, we have set upper limits on other decay modes of this state. 
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FIG. 5: The mass difference spectrum M(D s tt°j) — M(D s ir°) for candidates for the decay 
D s j(2463) — > D*j(2317)7, after application of the selection criteria described in the text. 



We have measured the mass splitting of this state with respect to the D s meson to be 
350.0 ± 1.2 [stat.] ± 1.0 [syst.] MeV/c 2 , and we find its natural width to be T < 7 MeV at 
90% C.L. 

We have observed and established the existence of a new narrow state with a mass 
near 2.46 GeV/c 2 in its decay to D* + tt°, which we have denoted D s j (2463). We have 
demonstrated that the signal for this decay cannot be interpreted as a reflection from the 
-D*j(2317) + — > D+n decay. The measured properties of this state are consistent with its 
interpretation as the 1 + partner of the + state in the spin multiplet with light quark angular 
momentum of j = 1/2. We have measured the mass splitting of this state with respect to 
the D* meson to be 351.2 ± 1.7 [stat.] ± 1.0 [syst.] MeV/c 2 . The natural width of this state 
is found to be T < 7 MeV at 90% C.L. Since the -D s j(2463) mass lies above the kinematic 
threshold for decay to DK (but not for D*K), the narrow width suggests this decay does 
not occur. Since angular momentum and parity conservation laws forbid a 1 + state from 
decaying to two pseudoscalars, this provides additional evidence for the compatibility of the 
D sJ (2A63) with the J p = 1+ hypothesis. 

In the model of Bardeen, Eichten and Hill [15], a J p = 1 + state is predicted with the 
same mass splitting AM with respect to the 1~ state as that between the + and 0~ states. 
Taking the difference between the two mean mass differences reported above, we obtain 
*(AM) = (351.2 ± 1.7) - (350.0 ± 1.2) = 1.2 ± 2.1 MeV/c 2 for the difference between the 
1 + — 1~ and + — 0~ mass splittings, where the dominant uncertainty is due to statistics. 
Thus our observations are consistent with these predictions. 
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